Delocalization in two dimensional disordered Bose systems and depinning transition 

in the vortex state in superconductors. 
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We investigate two-dimensional Bose system with the long range interactions in the presence 
of disorder. Formation of the bound states at strong impurity sites gives rise to an additional 
depletion of the superfluid density . We demonstrate the existence of the intermediate superfluid 
state where the condensate and localized bosons present simultaneously. We find that interactions 
suppress localization and that with the increase of the boson density the system experiences a sharp 
delocalization crossover into a state where all bosons are delocalized. We map our results onto the 
three dimensional system of vortices in type II superconductors in the presence of columnar defects; 
the intermediate superfluid state maps to an intermediate vortex liquid where vortex liquid neighbors 
pinned vortices. We predict the depinning transition within the vortex liquid and depinning induced 
vortex lattice/Bose glass melting. 

PACS numbers: 73.23Hk, 73.22Lp, 71.30.+h 



Effects of disorder on strongly correlated systems are 
central to contemporary condensed matter physics. Dis- 
ordered Bose systems are of special interest Q y, |^ . Su- 
perfluid '*He and '^He remain one of the favorite and im- 
portant experimental tools for studies macroscopic quan- 
tum effects. Moreover, quantum mechanical mapping 
that relates quantum Bose system to tfiermodynamics 
of superconducting vortices in the space of the one di- 
mension up, motivates research of Bose fluids in random 
environment by the hope for further progress in vortex 
physics. 

There has been a significant recent advance in un- 
derstanding disordered Bose systems started by 0, 0| 
who discussed a continuum model of the dilute inter- 
acting Bose gas in a random potential. The proposed 
model described the quasiparticle dissipation and deple- 
tion of superfluidity at zero temperature. Later we de- 
veloped a systematic diagrammatic perturbation theory 
for the dilute Bose gas with weak disorder that enabled 
us to extend the description of superfluid thermodynam- 
ics on finite temperatures Disorder corrections to the 
thermodynamic potential and the disorder-induced shift 
of the condensation temperature resulting from disorder 
scattering of quasiparticles were found. In particular, we 
obtained that the superfiuid density decreases monoton- 
ically with the temperature. 

Yet there are many questions that remain open. First, 
more close examination of models IS shows that 

the described depletion of superfluid density reflects only 
scattering of zero energy quasiparticles by random poten- 
tial, while possible contribution from the quasiparticles 
bound states was overlooked. Further, the fact that some 
quasiparticles leave condensate and get localized suggests 
the existence of the novel intermediate state where super- 
fluid and localized components present simultaneously. 
This poses the next question: if this intermediate state 
exists, what part of the phase diagram does it cover? In 



other words the question is: would even arbitrarily weak 
disorder localize the part of the condensate, or else this 
localization effect vanishes if disorder is too weak or if 
the Boson interactions are sufficiently strong? 

In this Letter we examine the behavior of 2D inter- 
acting Bose-system subject to strong disorder having in 
mind application of our results (via quantum mechanical 
mapping) to the system of 3D superconducting vortices 
in type II superconductors pinned by columnar defects. 
It is well established that at low temperatures where the 
pristine vortex structure would form Abrikosov lattice, 
columnar defects cause formation of strongly pinned Bose 
glass At the melting Hne Bm{T), Bose glass melts 

into a vortex liquid. In the related quantum mechani- 
cal 2D system vortex liquid corresponds to a superfluid 
Bosons. The depletion of the superfluid density by disor- 
der corresponds to enhancement of the average tilt modu- 
lus of the related vortex system 0,0] . Moreover, one can 
expect that this stiffening is accompanied by the change 
in transport properties: the part of the vortices remain 
pinned and does not participate in transport. This may 
be interpreted as the intermediate vortex phase where 
both vortex liquid and pinned vortices coexist. 

We focus on the most interesting case where the ap- 
plied magnetic filed is not too close to Hd such that dis- 
tance between the vortices is less than their interaction 
range, which coincides with the magnetic field penetra- 
tion length A. We show that upon increasing the temper- 
ature (or, equivalently, magnetic field ), the system un- 
dergoes a sharp crossover between the phase containing 
vortices pinned by the columnar defects and the phase 
where all vortices are delocalized and find the expression 
for the crossover line: 
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where p is the density of vortices, K = A/^^lis the ra- 
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tio of the magnetic penetration and coherence lengths, 
is the characteristic size of the columnar defect, e is the 
anisotropy parameter, the length It is determined by the 
temperature T as It = (^'o/47r)^/T, and c is a numeri- 
cal coefficient of the order one. The localization length 
L that enters the logarithmic factor is defined below in 
Eq. (|10|l . The corresponding phase diagram is shown on 
Fig. 1 where the delocalization line is drawn along with 
the vortex lattice melting line. In what follows we formu- 
late our model in terms of vortex configuration, establish 
the existence of the delocalization crossover, and derive 
equation Q for the crossover line. 

Model. A free energy of the vortex configuration can 
be written as a functional of the vortex displacements 
ri{z) as 



F 



dz 



where the interaction between the vortices is 



v{r)^eoKoir/X), 



(2) 



(3) 



with eo = ($o/47rA)^ and Kq being the Bessel function. 
The vortex tension ei is related to the anisotropy param- 
eter e and eo as ei = e^eo and the disorder potential U{r) 
can be written as a sum of the potential wells represent- 
ing the columnar defects 



U{r) = u{r - 



(4) 



where is the coordinate of the i — th defect. 

Statistical mechanics of vortices can be mapped onto 
the quantum mechanical problem of interacting two- 
dimensional bosons described by the Hamiltonian 

H = j d^r^'^{r)[p^/2m- n + U{r)]^{r) 

+ ^ j d^rid^r2n{ri)v{ri-r2)h{r2), (5) 

where ^\ are bose creation and annihilation operators, 
/i is the chemical potential and n{r) = ?/'t(r) V'(^) is the 
density operator of the bose gas. Results obtained in 
the Bose gas representation can be translated onto vor- 
tex language by the substitute fi T, m —^ ei. The 
energy of the Bose gas corresponds to the free energy per 
unit length of the vortex lattice. Keeping in mind the 
application to vortices we will make several assumptions 
about the form of a single potential well: We will consider 
temperatures T « Tc, in this regime pinning potential is 
given by 
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where ^ is the coherence length and rr is the scale 
that characterizes the size of the columnar defect. The 
strength of the pinning potential can be characterized by 
the dimensionless coefficient 



(7) 



We will assume that pinning is weak such that /3 <C 1. In 
this limit the solution of the Schrodinger equation 

[-Vr/2ei + u{r)]^{r) ^ El ijir), (8) 

results in the pinning energy 

rp2 
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The localization length of the quantum sate is related 
to the bound state energy as ^ /\Ei\m, translating 
this relation on vortex language we define the localization 
length of the vortex as on the columnar defect as 



(10) 



Further, we assume that the interaction between the 
bosons (vortices) is strong enough such that the bound 
state of the well can be only single occupied. This condi- 
tion is satisfied when the energy of the double occupied 
state E2 is much larger than the the absolute value of the 
energy of the single occupied state Ei . The interaction 
energy of the double occupied state is 

E2 = \ Jd^rid\2 ri{ri)Mri)v{ri-r2)riir2)Mr2), 

(11) 

where (f>i is the wave function of the bound sate. In the 
case when the localization length L is smaller than the 
interaction range A the energy E2 can be estimated as 



E2 ~ eo ln(A/i), 



(12) 



where we have used the asymptotic of the Bessel function 
Kq{x) w ln(l/a;). The condition of no double occupation 
can be now written as 



eoei ln(A/L) » {T/Lf 



(13) 



In case of low concentration of defects we can consider 
the contributions of different potential wells separately, 
thus we end up with the problem of the interacting Bose 
gas and single potential well ( or a single columnar defect 
placed in the vortex liquid). Leaving only one potential 
well we can easily diagonolize the non-interacting Hamil- 
tonian and formulate the problem in terms of eigenfunc- 
tions of the noninteracting Hamiltonian 



(14) 



where (j)k are the eigenfunctions of the non-interacting 
Hamiltonian 



[f/2m + u{r)] (t>k{r) = Ek M^), 



(15) 
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and Ek is the energy of the k — th state. The bound state 
corresponds to the index k = 1. Although this state has 
the lowest energy, the Bose-condensation to this state 
cannot occur due to our assumption of no-double occu- 
pation. Thus, the condensation takes place to the ex- 
tended state with the lowest energy (which we label with 
k — 0). The normalization conditions for eigcnfunctions 
(j)k are as follows: 



'k9k = 



V, 
1, 



Ek >-0 
k = 1 



(16) 



with V being the volume. In the absence of the current 
the basis can be chosen real, thus we will assume further 
on that (p* = (p. The '0-operators that enter the Hamil- 
tonian (jsj can be presented as a sum of three terms 



(pk bk 



(17) 



Sfc>0 



representing the condensate, the bound state, and all 
other states having energy larger than the condensate 
energy. The effective temperature of the Bose gas is 
related to the inverse physical length of the supercon- 
ducting sample in the z direction. We will assume that 
the size of the sample is microscopically large, this cor- 
responds to zero temperature limit for bosons. We will 
consider the most interesting case where the applied mag- 
netic filed is not to close to Hd, such that the average 
distance between the vortices is smaller than the screen- 
ing length A. In this regime the mean filed Bogolubov 
approximation, that we will use, is justified for the case 
of sufficiently high density of vortices. [l3| Thus, the 
third term in Ea. (|17(l representing the particles that do 
not belong to the condensate is small and can be omitted 
in the leading order. 

Now we turn to the analysis of the population of the 
bound state in the presence of the condensate. Inserting 
the representation Ea. (|17|l into the Hamiltonian Q and 
keeping only two first terms in the representation (|17|l we 
get the effective Hamiltonian 



H, 



b\{Ei+a~ ^,)bl+J{b\+h), (18) 



where the terms b\b[ were not written due to no double 
occupancy condition. The coefhcients a and 7 in the 
Hamiltonian H18|l are 

a = bl J d\id\2[(pl{ri)v{n~r2)<j)l{r2) 

+ (t>i{ri)(l)o{ri) v{ri - r2)(f)i{r2)(f>o{r2)], (19) 
7 = 6^ / d\id\2cl>Un)v{ri-r2)cl>{r2), (20) 

where the average value bo of the operator 60 is related 
to the condensate density po as ~ po- The model (fH^ 
can be easily solved: Presenting the wave functions in 
the form 



iP = ao\0 > +ai\l >, 



(21) 



with oq and oi being the amplitudes of the zero and single 
occupied states we find two eigenstates corresponding to 
the energy levels 



E± = 



E ± y^E^ + 472 



(22) 



where E = Ei + a — fi. We see that E+ > and E^ < 
for any sign of the energy _B, thus the state with energy 
is always occupied while E^ is always empty. The 
occupation of the center is determined by the state with 
the lowest energy (E^) and is given by 



n = fli = 



{EM^ + {Eh),/{E/W 



(23) 



One can easily see that n — s- 1 when E/j <ti —1 and n 
when iJ/7 ^ 1. Thus the quantity E/'^ controls the oc- 
cupation such that a,t E = localization-delocalization 
crossover takes place. Making use of Ea. (|19|l and rela- 
tionship between the chemical potential, condensate den- 
sity and interaction p — 6§uo where vq — J d^rv(r), for 
the parameter E we find E = Ei + SE with 



SE 



d^ri d^r2 0i(ri) v{ri-r2) cj)i{r2) ~ pL^ co ln(A/L) 

(24) 



Now, using the relationship between the bound state 
energy and localization length, from the condition 
E = 0, we get the vortex density at the localization- 
delocalization crossover 



P = 



(25) 



where c ~ 1 is a numerical constant. Plugging in the 
localization length from Eq. H1U|I we arrive at the expres- 
sion for the delocalization crossover line which can be 
rewritten in terms of the depinning field for vortices as 



B, 



dep 



$0 



27rcr2 



eoei^HX/L) 



exp 



T 



(26) 



where To is the effective temperature dependent depin- 
ning energy Tq = i/cieo/S. One can easily check that 
the assumed condition of no double occupancy (|13() al- 
ways satisfies near the crossover as long as pi? < 1. 

As we see from Ea. (|23f) . the occupation of the center 
changes continuously as long as "hybridization parame- 
ter" 7 is larger than zero. Therefore the phase boundary 
between the localized and delocalized phases should be 
viewed as a crossover rather than the true transition line. 
Note, however, that in case of finite concentration of de- 
fects the delocalization transition may result in signifi- 
cant change of the condensate density, that may induce 
the true superfluid-bose glass transition studied in 0| . 

While the ultimate type of the crossover remains an 
open question, experimentally it can appear very sharp 
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Vortex liquid coexisting i 
with pinned vortices i 




Delocalized vortex liquid 



Delocallzation crossover line 



Expected melting line 



FIG. 1: Schematic phase diagram for the vortex system with 
columnar defects. The melting line under assumption of pin- 
ning is denoted by Bm, while the melting line of the pristine 



(0) 



lattice is denoted by B 



and undistinguishable from the true first order transition. 
Indeed, the sharpness of the crossover is controlled by the 
ratio of 7 and 6E. Formula (|20l) determining 7, should be 
corrected to include the screening of the potential by the 
surrounding Bose gas. In the dense limit the screening 
length is given by Igc ~ A(l/27rpe^ l^y^^. The parameter 



7 is estimated then as 



7 ; 



pLX^eo 



It£ V 27r 



and the ratio "f /SE becomes: 

-f/SE - AV(/Tie ln(A/L)) 



(27) 



(28) 



and is always small due to the no-double occupancy con- 
dition We thus conclude that the crossover between 
the localized and non-localized phases is sharp. 

We are now in a position to construct a phase diagram 
for the disordered Bose system. We use the correspond- 
ing 'superconductor vortex language'; the schematic vor- 
tex phase diagram for the system with columnar defects 
is presented in Fig 1. Shown are the melting line of the 
pristine sample B^^ , and melting line Bm in the sample 
with the columnar defects shifted, as it should 0,01 '^P" 
wards as compared to . Note that because of the ex- 
ponential dependence of temperature, Bdep drops much 
faster than Bm as temperature approaches Tc- Thus the 
depinning line that lies in the vortex liquid domain at 
high sufficiently fields hits and then crosses the melting 
line Bm , traversing then to Tc and crossing pristine melt- 
ing line at very low fields. 

In the liquid phase the depinning line marks the 
crossover separating the intermediate state with the 
columnar defects occupied by vortices (i.e. vortex liq- 
uid coexists with the pinned vortices), and the phase 



where all vortices are delocalized from the columnar de- 
fects. In the strip confined between the Bm'' and Bm 
lines, Bdep marks melting of the vortex lattice. Indeed, at 

si? < B < Bm lattice is stabilized by vortices pinned by 
columnar defects. As soon as depinning at Bdep occurs, 
columnar defects become inessential, and lattice looses 
its stability since in the absence of defects vortex liquid 
is a stable thermodynamic state above Bm^. We there- 
fore discover a novel type of the vortex lattice melting, 
the depinning induced melting, that occurs in the interval 
Bm^ < B < Bm- 

This new kind of vortex lattice melting may well ex- 
plain the origin of the low-field kink in the melting line 
that clearly indicates a switch between the different melt- 
ing mechanisms. This kink has been observed in almost 
all experiments on the Bose glass melting and was 
addressed specifically in the recent study of the melting 
of "porous" vortex matter Although these results 

obtained for pronounced vortex lines may be applied only 
with caution and many reservations to BSCCO, it is in- 
teresting to note that at temperatures around 80K the 
characteristic energy Tq for BSCCO parameters can be 
estimated as Tq ~ 5K, close enough to that observed in 

m- 

In conclusion we have shown that in a disordered Bose 
system with the long range interactions an intermedi- 
ate phase where both superfluid and localized bosons 
can exist. We demonstrated that a sharp delocallzation 
crossover occurs with the increase of bosons density. We 
apply our results to description of a vortex system in type 
II superconductors in the presence of columnar defects 
and show that this delocallzation crossover describes the 
depinning line separating two liquid phases, the interme- 
diate phase containing pinned vortices and 'fully molten' 
liquid. We predict a new kind of the first order vortex 
lattice melting, depinning induced melting, at low fields. 
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